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We investigate the synthesis of continuous-variable two-mode unitary gates in the setting where 
two modes A and B are coupled by a fixed quadratic Hamiltonian H . The gate synthesis consists of 
a sequence of evolutions governed by Hamiltonian H interspaced by local phase shifts applied to A 
and B. We concentrate on protocols that require the minimum necessary number of steps and we 
show how to implement the beam splitter and the two-mode squeezer in just three steps. Particular 
attention is paid to the Hamiltonian xaPb that describes the effective off- resonant interaction of 
light with the collective atomic spin. 



PACS numbers: 03.67.-a, 42.50.Dv 



I. INTRODUCTION 

One of the central problems of the quantum informa- 
tion theory is to establish what resources are sufficient 
for universal quantum computation. In this context, the 
question whether a given Hamiltonian H can simulate 
another one has attracted considerable attention recently 
[HQ EI [J 00 In its simplest form, this problem may 
be formulated as follows. Consider two parties, tradition- 
ally referred to as Alice and Bob, possessing a single qubit 
each. The interaction between those two qubits is gov- 
erned by a fixed Hamiltonian H, that is determined by 
the physical properties of the systems that represent the 
qubits. In addition to the interaction H , Alice and Bob 
may attach (local) ancillas to their qubits and perform 
arbitrary local unitary operations on their subsystems. 
It is usually assumed that these local operations are very 
fast compared to the evolution induced by the Hamilto- 
nian H. The task for Alice and Bob is to simulate the 
evolution due to different Hamiltonian H' . Two kinds 
of simulations should be distinguished. The infinitesimal 
time simulation Q, d H, IE El consists of simulating the 
action of the Hamiltonian H' for an infinitesimally short 
time interval At. The gate synthesis 0, H ES 111 113 
requires the implementation of the unitary transforma- 
tion U' — exp(—iH't) for some finite time t. 

It turns out that in the two-party setting all nonlo- 
cal Hamiltonians H are qualitatively equivalent. Given 
enough time r, Alice and Bob can, with the help of local 
ancillas, simulate the evolution exp(-iH't) for any H' 
J2). The central question, then, is what is the optimal 
simulation. The latter may be defined as a simulation 
that requires the shortest time. For the two-qubit case, 
this problem has been completely solved and the optimal 
protocols for Hamiltonian 0, @ and unitary gate [g sim- 
ulations have been determined. The situation becomes 
much more complicated for higher dimensional systems 
and for higher number of involved parties. Simulation 
protocols suggested for these generic settings are rather 
involved and it is not known which protocols are optimal. 

It should be stressed that most of the work focused on 
discrete variable systems: qubits or, more generally, qu- 



dits. Recently, however, Kraus et al. extended the notion 
of Hamiltonian simulation to continuous variable systems 
|l3j. They assumed that Alice and Bob possess a single- 
mode system each and these two modes are coupled via 
quadratic Hamiltonian (we assume Ti = I throughout this 
paper): 

H = cux A x B + ci 2 x A p B + C21PAXB + C22PAPB, (1) 

where Xj and pj are two conjugate quadratures of the jth 
mode. Kraus et al. showed that almost every Hamilto- 
nian Q is capable to simulate any other Hamiltonian of 
the form Q provided that Alice and Bob can apply fast 
local phase shifts described by single-mode Hamiltonians 



Ha = x 2 A + p A and H B 



+ Pb- 



These results are interesting both from the theoretical 
and experimental points of view. In particular, the off- 
resonant interaction of ligh t with the collective atomic 
spin [3 [H [H IHIli [la E3 can be described by the 
effective unitary transformation 



U = exp(—itH A L), 



where the Hamiltonian 



Hal 



kxaPb 



(2) 



(3) 



is a special instance of (Q. The typical geometry of the 
experiments is such that a light with strong coherent field 
polarized along the x axis propagates along the z axis 
through the atomic sample, whose spin is also polarized 
along the x axis. The x and p quadratures are defined 
as the properly normalized y- and z-components of the 
collective spin operators describing the polarization state 
of light and atomic ensemble, respectively 0, [l6L Il7| . 
In recent beautiful experiments it was demonstrated that 
the interaction J2J can be employed to squeeze the atomic 
spin E3 , entangle two distant atomic ensembles E3 1 an d 
transfer the quantum state of light into the atoms E3- 
Schemes for teleportation and swapping of the quantum 
state of collective atomic spin have been suggested 0, 
Il7j . These experiments and proposals in fact rely on 
the quantum non-demolition (QND) measurement of the 
atomic quadrature, possibly accompanied by a suitable 
feedback. 
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As showed by Kraus et al, the Hamiltonian © can 
simulate any Hamiltonian Q. In particular, Hal can 
be used to implement a beam splitter and a two-mode 
squeezer. This is very appealing since it suggests that, 
for instance, the storage of the quantum state of light into 
atoms and the subsequent readout of the quantum mem- 
ory - the transfer of quantum state of atoms onto light - 
can be implemented in a unitary way if the Hamiltonian 
(|3"|) is used to simulate a beam splitter. 

However, currently there are technical difficulties that 
will complicate the actual practical realization of this 
procedure. The effective unitary transformation J2J de- 
scribes the modification of the polarization state of the 
light pulse after the passage through the atomic sample 
due to its coupling with atoms. This means that the 
Hamiltonian simulation requires several passages of the 
light pulse through the atomic sample (c.f. the detailed 
description of the simulation protocol in Sec. II). In cur- 
rently envisaged experiments, the pulse width must be 
at least 1 /zs |2l| which corresponds to the length 300 
m. This long pulse would have to be stored somewhere 
(e.g. in an optical fiber) until its tail leaves the atomic 
sample. Only then can the pulse (with properly applied 
phase shifts) be fed to the atomic sample again. 

These practical considerations imply that the approach 
relying on the infinitesimal time simulation is not very 
convenient from the experimental point of view. It is 
possible to simulate a gate by concatenating a large se- 
quence of short-time Hamiltonian simulations but this 
would require a large number of manipulations and pas- 
sages of the light pulse through the sample. Since, in 
practice, every round of the gate synthesis procedure is 
necessarily accompanied by some losses and other errors, 
the accumulation of the errors would negatively influence 
the simulation. 

In this paper, we show how to implement several im- 
portant two-mode interactions with the Hamiltonian 
such that the number of the applications of the Hamilto- 
nian JQ) is minimized. We demonstrate that only three 
sequences of evolution governed by Hamiltonian H , inter- 
spaced by (fast) local phase shifts on both subsystems, 
suffice to implement a two-mode squeezer and a beam 
splitter. For the specific Hamiltonian © we also design 
a single-mode squeezing gate that involves four evolu- 
tion steps. These results illustrate that several important 
quantum information processing tasks such as entangling 
the light and collective atomic spin, or a transfer of the 
quantum state of light into atomic clouds and vice versa, 
can be carried out with a small number of repeated pas- 
sages of the light pulse through the atomic sample. 

This paper is structured as follows. In Sec. II we in- 
troduce the notation, the canonical form of the interac- 
tion Hamiltonian ^ and we describe the gate simulation 
protocol. In Sec. Ill, we consider the simple interaction 
Hamiltonian and we show how to implement the two- 
mode squeezing operation, beam splitter transformation 
and also single-mode squeezing as a sequence of three (or 
four) intervals of evolution governed by Hamiltonian (J3J) 



combined with local phase shift operations. In Sec. IV 
we extend this analysis to the generic interaction Hamil- 
tonians Finally, Sec. V contains conclusions. 

II. DESCRIPTION OF THE SIMULATION 
PROTOCOL 

In this section we describe the simulation protocol. 
The gate synthesis consists of a sequence of N intervals 
of evolution governed by the Hamiltonian H followed by 
local unitary phase shift transformations. The resulting 
unitary gate G is given by 

G = Vle- lHtN V N . . . Vje~ iHt2 V 2 V?e~ lHtl V 1 . (4) 

The local phase shift operation applied to modes A and 
B reads 

Vj = e -' l ^ Aaia g! e~ i ^ Bbib , (5) 

where a and b are the annihilation operators of modes A 
and B, respectively. Note that Vj +1 Vj is still of the form 
©, with <f) A = 4>.j A - <t>j+i,A and <j) B = (pjB - <t>j+i,B- 
With the help of the useful identity 

exp(—iHt)U = exp (-iU^HUt) (6) 

we can rewrite Eq. Q as 

Q _ p-i-tfivtjv e ~iH2t2 e ~ iHiti 

where Hj = VjHVj. 

The Hamiltonian is characterized by four parame- 
ters. However, by means of local rotations, we can always 
transform this Hamiltonian to a diagonal form 

H c = cix A p B + c 2 p A x B , (8) 

where c\ = a\ and c 2 = <j2<iet[C]/\ det[C]|, and o\ and a 2 
are the singular values of the matrix C defined as (C)y = 
Cij 0|- I n close analogy to the qubit case 0], we may 
refer to H c as the canonical form of H. Mathematically, 
we have 

cxp(-iH c t) = y f cxp{-iHt)V, (9) 

where V is a local rotation 101. This shows that, without 
loss of generality, we may assume that H has the canon- 
ical form (jHJ. In particular, it follows that H is able to 
simulate an arbitrary H' Q if and only if H c is able to 
simulate an arbitrary canonical Hamiltonian 

In Eq. (jSJl the phase shifts <f>j may be arbitrary. In 
what follows, we focus on the phase shifts that preserve 
the canonical form of H. There are four inequivalent 
possibilities: 
(a) (\> A = 0, <f) B = 0, 

Hi = cixaPb + c 2 p A x B . (10) 
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(b) <j} A = 7T/2, 4> B = 37T/2, 

iJ 2 = C.2XAVB + CiPAXB- (11) 

(c) ^4 =tt, (j) B = 0, 

i?3 = -CiZAPB - C2PAXB ■ (12) 

(d) A = tt/2, 0b - it/2, 

Ha = -c 2 xapb - CiPAXB- (13) 

From the structure of these Hamiltonians we can deduce 
that two different noncommuting canonical Hamiltonians 
Hi and H 2 are available. Furthermore, we can see that 
H3 = —Hi and H4 = —H2, hence we can implement any 
transformations of the form exp(— iHit) and exp(— iH 2 t) 
where t is an arbitrary real number, positive or negative. 
The two specific cases ci = c 2 and c\ — — c 2 when Hi = 
±i?2 and the simulation is not possible correspond to 
the Hamiltonians of a two-mode squeezer and a beam 
splitter, respectively. 



Without loss of generality, we may assume that the cou- 
pling constant is equal to unity, hence the two relevant 
Hamiltonians read 

Hi = x a Pb, H 2 = PaXb- (15) 

All the canonical Hamiltonians JHJ have the important 
property that the x and p quadratures are not mutually 
coupled when we write down the Heisenberg equations of 
motion for Xj and Pj. This means that the evolution of 
the operators x — (xa,xb) T and p — {pa,Pb) T is gov- 
erned by the following linear canonical transformations: 

x out = Sx in , p ou t = Rpin- (16) 

This decoupling of x and p quadratures greatly simplifies 
the analysis. The transformation l|16(l must preserve the 
canonical commutation relations [xj, pf.] — iSjk- From 
these conditions we can express the matrix R in terms of 
S, 

R=(S~Y, (17) 



III. XP COUPLING 

Having established the notation and described the sim- 
ulation protocol, we may proceed to the unitary gate 
synthesis. Namely, we would like to decompose the uni- 
tary transformation G that we want to simulate into a 
sequence of unitary evolutions governed by Hamiltonians 
Hi&nd H2 that were defined in the previous section, 

Q _ e -iH 2 t N e —iHitN-i e -iH 2 t 2 e ~iH 1 t 1 

We are particularly interested in the simulations that in- 
volve the lowest possible number of steps N, because such 
simulations require low number of local manipulations in 
the eventual experimental implementation. 

We note here that Eq. I|14|) is an example of a de- 
composition of a group element into a product of N 
other group elements. In the present case, the under- 
lying group is the symplectic group Sp(4, R) of all linear 
canonical transformations of the quadratures of the two 
modes A and B [22I I23I ] . It is worth mentioning here that 
the related problem of a decomposition of the symplec- 
tic transformation into a sequence of simple evolutions 
associated with the common passive and active linear 
optical elements has been studied recently. Braunstein 
has shown that any N mode symplectic transformation 
can be implemented as a sequence of an iV-mode passive 
linear interferometer followed by N single-mode squeez- 
ers and another passive interferometer — the so-called 
Bloch- Messiah decomposition j^. The decompositions 
of this kind have also been applied to investigate the 
properties of nonlinear optical couplers [2^, |2(j . 

In this section, we shall consider the simplest and also 
the experimentally relevant coupling between the two 
systems described by the interaction Hamiltonian JJJl. 



hence the evolution of p quadratures is uniquely deter- 
mined by the evolution of the x quadratures. 

Our task is to implement two-mode unitary gates 
(symplectic transformations) as a sequence of a small 
number of unitary transformations generated by the 
Hamiltonians (|15fl . The matrices Si and S2 associ- 
ated with the unitary evolutions Ui = exp(-iHit) and 
U2 = exp(—iH 2 t) read 

SiW = (j?), s rt) = ( 1 oi)- ( 18 ) 

The factorization Ijl4(l can be rewritten in terms of the 
matrices Sj as follows, 

S = S 2 (t N )Si(t N ^i) . . . S 2 (t 2 )Si(ti), (19) 

where S is the matrix associated with the gate G. Since 
det Si = dct f>2 = 1, we are restricted to a three paramet- 
ric subgroup of transformations S such that det S = 1 . 
In what follows, we will discuss the implementation of 
three important gates: a beam splitter, a two-mode 
squeezer and a single-mode squeezer. 



A. Beam splitter 

The beam splitter operation is described by the matrix 

S BS (ct>) = ( C ° S t Sin lV (20) 
y— sin0 cos0 J y ' 

We show that this transformation can be implemented 
as a sequence of three evolutions (1181) . 

SbsW ^ Si( 1 )S 2 (fJ)Si(a). (21) 
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The explicit multiplication yields 

8 »°M=( a +\Vtaft) !+%)• (22) 

If we compare the elements of the matrices on left- and 
right-hand sides of Eq. \22\. we obtain a set of equations 
for the parameters a, (3 and 7 whose solution yields: 

a = — tan— , p = sin0, 7 = a. (23) 

The parameters a and (3 can be negative but this is not 
an obstacle as explained in the previous section since we 
can change the sign of the Hamiltonian Hi or H 2 by 7r 
rotation of one of the systems. Two cases of particular 
importance are (i) the balanced beam splitter (0 = 7r/4), 
that requires a = 1 — y/2 and (3 = y/2/2, and (ii) the 
swap (4> = 7r/2) that exchanges the quantum states of 
the two systems, a = — 1, /3 = 1. 

The swap gate is closely related to the two-step proto- 
col for mapping the state of collective atomic spin on light 
that was suggested by Kuzmich and Polzik [20]. In fact, 
their two-step protocol can be obtained by simply remov- 
ing the last step of the present three-step swap gate. The 
removal of the third step means that the mapping adds 
some noise and the procedure is thus only approximate. 
A possible way how to improve its performance is to use 
squeezed light. For details, see [2Cj . 

B. Two- mode squeezer 

Let us now turn our attention to the two-mode 
squeezer, described by the following matrix, 

„ , , ( coshr sinhr \ 
5 ™ s W=( v sinhr cosh. J' ^ 

Similarly as in the case of the beam splitter, we attempt 
to implement this transformation as a sequence of three 
evolutions, c.f. Eq. i|21|) . By comparison of the right- 
hand side of Eq. I)22|l with the matrix l|24|l we again 
obtain a system of nonlinear equations for the parameters 
a, f3, and 7 having the solution 
v 

a = tanh — , (3 = sinhr, 7 = a. (25) 

Note the similarity with the results for the beam splitter, 
the difference is essentially that the goniomctric functions 
have been replaced by their hyperbolic counterparts. The 
parameters are finite for any finite r. However, /3 grows 
exponentially with r and for large r we have (3 cx e r . On 
the other hand, for small r we get (3 f=a r. This implies 
that we may reduce the synthesis time if we implement 
the two-mode squeezing transformation as a sequence of 
n two- mode squeezers with r' — r/n. The reduction of 
the time is achieved at the expense of a higher number of 
steps of the gate synthesis protocol. For modest values 
of squeezing r, the sequence of three evolutions is advan- 
tageous since it involves the minimum necessary number 
of manipulations of the systems A and B. 



C. Single-mode squeezer 

After dealing with two-mode gates, let us now focus on 
the single-mode gate, namely the single-mode squeezer. 
Since det S = 1 , the squeezing of quadrature xa will 
necessarily by accompanied by anti-squeezing of xb , and 
vice versa, 

S SMS (r)= (0 e ° r ). (26) 

It turns out that in this case a sequence of three trans- 
formations l|18fl is insufficient and we must consider a 
sequence of four basic evolutions: 

SsmsM = S 2 (S)S 1 ( 1 )S 2 (/3)S 1 (a). (27) 

We proceed as before and derive equations for the four 
parameters appearing in (|27ll : 

e r = 1 + a[3 + 6 (a + 7 + a(3j), 

e- r = 1 + 7/3, 
= f3 + 5(l + 1 (3), 
= a + 7(1 + a/3). 

This system of equations has a one-parametric class of 
solutions, given by 

e r - 1 e r (l - e r ) 

13 =- 1 = -ae~\ 6= e[L & \ (28) 

a a 

and a is arbitrary but nonzero. We may choose the op- 
timal value of a that minimizes the total time T needed 
for the implementation of the squeezing operation, 

T = M + + | 7 | + |*|. (29) 

Assuming that r > we obtain by solving dT/da = 
the optimal value, 

^V— < 30 > 

In the limit of small r, all four parameters a, f3, 7, and 
5 are proportional to y/r. This stems from the fact that 
the single-mode squeezing Hamiltonian H$ms = xaPa 
cannot be obtained as a linear combination of the two- 
mode Hamiltonians H\ and H 2 and only the terms of 
the order of 0(t 2 ) or higher in l|14l) may give rise to the 
contribution proportional to H$ms- 

IV. GENERIC QUADRATIC COUPLING 

In this section we shall assume that the interaction 
Hamiltonian has the generic canonical form. Although 
the mathematical analysis will be more involved we shall 
still be able to derive analytical formulas for the inter- 
action times characterizing the gate synthesis. Without 
loss of generality, we may assume that c\ = 1 in |JHJ|. We 
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have to distinguish two classes of Hamiltonians giving rise 
to qualitatively different evolutions of the quadratures in 
the Heisenberg picture. For c 2 > the dynamics resem- 
bles an amplifier while for c 2 < we obtain oscillatory 
dynamics reminiscent that of a beam splitter. We shall 
discuss these two cases separately. 



A. Amplifier-like Hamiltonians 

Suppose first that c 2 > and introduce a more conve- 
nient notation c 2 = s 2 , s > 0, hence 

Hi = xaVb + s 2 p A x B , H 2 = s 2 x A p B + p A x B - (31) 

It is an easy exercise to derive matrices Si and S 2 corre- 
sponding to the unitary evolutions governed by Hi and 
H 2 , respectively, 



S+(t) 



cosh(si) s sinh(st) 

| sinh(si) cosh(st) 

cosh(st) - sinh(si) 

s sinh(st) cosh(si) 



(32) 



In what follows we will focus on the implementation 
of the beam splitter and two-mode squeezer transforma- 
tions. We have seen in the previous section that these 
transformations could be implemented as a sequence of 
three basic evolutions governed by Hamiltonians Hi or 
H 2 that were of the form l|15|). Moreover, there was an 
inherent symmetry in this gate synthesis; we have found 
that 7 = a. It turns out that these basic symmetry 
properties remain valid also for the generic Hamiltonians 
(|31[) . and we can thus decompose the two-mode squeezing 
transformation as 

S TM s(r) = S+(a/s)S+((3/s)S+(a/s). (34) 

Here the parameters a — sti and (3 = st 2 are the rescaled 
interaction times. The nonlinear equations for the pa- 
rameters a and (3 are much more complicated than be- 
fore. Nevertheless, analytical results can be obtained. 
From the condition S12 = S 2 i we get 



tanh/3 



2 sinh a cosh a 



cosh a — (s 2 + 1 + s 2 ) sinh a 



(35) 



The condition Sn = S 22 is satisfied due to the symmetry 
(7 = a), and the parameter a can be determined from 
the last independent equation Si 2 /Sn = tanh7', which 
yields 



2y(s + s- 1 ) 



i + ?/(.s- 2 + i + s 2 ; 



= tanh 1 



where y — tanh a. This is a quadratic equation for y 
whose solution reads 



tanh a 



+ s -1 - a/1 + (s -2 + 1 + s 2 )(coshr)" 
(s~ 2 + 1 + s 2 ) tanhr 



(36) 



We have selected the root that yields the correct limit 
tanh a = when r — > 0. In the opposite limit r — > 00 we 
obtain 



tanh a c 



1 



s- 1 + 1 + s 

On inserting this back into Eq. (l.> r >[) we find that 
lim tanh (3=1. 



(37) 



(38) 



It is easy to check that the equations for a and (3 have 
finite solutions for any finite r. In the limit r — > 00, a 
approaches a finite asymptotic value, cf. Eq. I|37|l . while 
(3 grows to infinity. 

Suppose now that we want to implement the beam 
splitter transformation (|20|l . The calculations of the 
parameters a and (3 parallel those for the two-mode 
squeezer. From the condition Si 2 = — S21 we express 
[3 in terms of a: 



tanh/3 



-2 sinh a cosh a 



cosh 2 a 



j 2 — 1 + s 2 ) sinh a 



(39) 



(33) Since (s~ 2 - 1 + s 2 ) > 1, it follows that |tanh/3| < 



|tanh(2a)|. From the condition S12/S11 = tan< 
tain quadratic equation for tanh a leading to, 



we ob- 



tanh a 



(s 1 - s) - sign(s 1 -s)^J~- 



(s 2 — 1 + s 2 ) tan (f> 



The sign function in the above formula selects the root 
that yields the correct limit a = when <p — > 0. The 
formula l|4(JI) is applicable in the interval </> £ [0,7r/2]. In 
the limiting case 4> = 7r /2 we have 



tanh a n /2 



sign(s — s 
Vs~ 2 - 1 4 



c2 ' 



(41) 



which implies that | tanh a v / 2 \ < 1 iff s ^ 1. Further- 
more, it can be shown that a is a monotonic function 
of 4> in the interval [0, n/2]. We can thus conclude that 
with the interaction Hamiltonian of the amplifier type 
(|31f) we can implement any beam splitter transforma- 
tion 1)20(1 with the mixing angles in the interval [0,7r/2], 
which includes the two important cases of a balanced 
beam splitter (<fi = ir/4) and the swap (4> — ir/2). It fol- 
lows from the expressions l(39JI and (|40|l that the simula- 
tion becomes more and more time consuming for Hamil- 
tonians close to the two-mode squeezing Hamiltonian 
Htms = xaPb +Paxb, i-e., when s —> 1. 



B. Beam splitter-like Hamiltonians 

Having derived the gate synthesis parameters for the 
Hamiltonians l|31|l . we proceed to the interaction Hamil- 
tonians leading to oscillatory dynamics, 



Hi = x A p B -s paXb, 



Ho 



-s 2 xaPb+Paxb- (42) 
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The S matrices associated with these Hamiltonians read 



cos(st) — ssin(si) 
- sin(st) cos(si) 



cos(st) - sin(si) 
-ssm(st) cos(si) 



(43) 



(44) 



We shall not repeat the details of the derivations of the 
parameters a and (3 and we only summarize the results 
here. The beam splitter transformation can be accom- 
plished by the following choice: 



tan/3 



-2 sin a cos a 



cos 2 a + (s 2 + l + s 2 )sin a 



(45) 



and 



tana = ± 8 - ± 1 + 1 , (46 ) 

(s- 2 + l + s 2 )tan0 

This reveals that simulation of any beam splitter with 
4> € [0, 7r/2] is possible and the parameters satisfy \a\ < 
tt/2 and \(3\ < tt/2. 

Consider now the two- mode squeezing operation (|24|l . 
After some algebra, one obtains 



tan/3 = 



2 sin a cos a 



cos 2 a — (s 2 — 1 + s 2 ) sin a 



(47) 



and 



tana 



— s — sign(s 1 — s) 



— I 

cosh r 



(s 2 — 1 + s 2 ) tanhr 



(48) 



The function tana must be real which implies that 
the term under the square root must be non-negative. 
This constraint, in turn, limits the amount of two- mode 
squeezing that can be produced via a three-step protocol 
(|34|l . It holds that r < r t h where 



coshr t h = y s 2 — 1 + s 2 



(49) 



The origin of this bound lies in the fact that the dynam- 
ics governed by the Hamiltonians i|42|) and captured by 
the matrices l|43[l and (|44|1 is oscillatory and fully peri- 
odic with period 2ir/s. Squeezing above the threshold rth 



can be achieved only if we concatenate several three-step 
protocols. It thus appears that the amplifier- like Hamil- 
tonians H31JI are in certain sense more versatile than the 
beam-splitter like Hamiltonians (|42|l . because the former 
allow to implement any two-mode squeezing gate and 
also any beam splitter with <f> G [0,7r/2] via a three-step 
protocol 



V. CONCLUSIONS 

In this paper we have addressed the problem of gate 
synthesis for continuous variable systems. We have as- 
sumed that two single-mode systems A and B interact 
via a quadratic Hamiltonian. We have studied how to 
implement a unitary symplectic gate G with the use of 
the interaction Hamiltonian H as a resource. The gate 
synthesis protocol consists of a sequence of evolutions 
governed by H and followed by fast local phase shifts 
applied to the systems A and B. We have focused on 
the gate simulation protocols that involve the minimal 
necessary number of steps, because these protocols re- 
quire a low number of local control operations which is 
important from the experimental point of view. We have 
shown that a three-step protocol suffices for simulation 
of the two-mode squeezer as well as a beam splitter. For 
the specific case of the Hamiltonian J3| we have also es- 
tablished a four step implementation of a single-mode 
squeezer. Our results are applicable to any physical sys- 
tems coupled via quadratic Hamiltonians. In particular, 
the gate synthesis protocols proposed in the present pa- 
per may find applications in the experiments where light 
interacts with atomic ensembles via a Kerr-like coupling 

Guam. 
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